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ABSTRACT 

Outflow angles and bulk Lorentz factors for 43 sources that have proper motions 
compiled by Vermeulen & Cohen (1994) are computed on the basis of Doppler 
factors and observed apparent motions in the plane of the sky. These estimates of 
outflow angles and bulk Lorentz factors are discussed along with their agreement with 
orientation unified models of AGN. 

Intrinsic (i.e. rest frame) brightness temperatures computed by using the inverse 
Compton and equipartition Doppler factors are discussed along with their relevance to 
the "Inverse Compton catastrophe" . Intrinsic luminosity densities and luminosities are 
presented, and the role of systematic errors is discussed. 

These studies are carried out using a sample of 100 compact radio sources compiled 
by Ghisellini et al. (1993). Error estimates for previously computed inverse Compton 
Doppler factors and equipartition Doppler factors are presented for these sources, 
along with a few updates of these Doppler factor estimates. 



Subject headings: BL Lacertae objects: general — galaxies: active — galaxies: 
kinematics and dynamics — quasars: general — radiation mechanisms: nonthermal 
relativity 



1. INTRODUCTION 



It is generally accepted that relativistic outflows from the core regions of radio-loud active 
galactic nuclei (AGN) are responsible for many of the interesting phenomena that are observed 
in these sources. Relativistic motion of synchrotron emitting plasma will result in the Doppler 
boosting of the synchrotron radiation from these outflows (discussed by many authors, e.g. 
Marscher 1987). In addition, relativistic outflows would readily explain the apparent superluminal 
motion observed on VLBI (Very Long Baseline Interferometry) scales (e.g. Pearson &: Zensus 
1987). 
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The Doppler factor is a key quantity in deriving intrinsic properties of an AGN. It is well 
known that the unbeamed flux density of a radiating source is related to the observed flux density 
by the Doppler factor. When combined with the apparent motion in the plane of the sky, the 
Doppler factor can be used to estimate the Lorentz factor and the viewing angle of an outflow, 
assuming the motion responsible for the Doppler boosting has the same speed as the motion of 
radio features in the plane of the sky (Ghisellini et al. 1993; Daly, Guerra, & Giiijosa 1997). The 
values of these Lorentz factors and viewing angles are significant in understanding AGN physics 
and are relevant to orientation unified models of AGN. 

There are several methods for computing the Doppler factor. One method, called the inverse 
Compton Doppler factor, is derived by assuming the observed X-ray flux is caused by inverse 
Compton scattering of synchrotron photons off the radiating particles (e.g. Marscher 1987). A 
second method, called the equipartition Doppler factor, is derived assuming that the sources are 
at or near equipartition of energy between radiating particles and magnetic field (Readhead 1994). 
Other methods (not discussed further here) include those based on variability measurements (Rees 
1967) or jet/counter jet flux ratios (Conway 1982). Ghiselhni et al. (1993; hereafter GPCM93) 
compute the inverse Compton Doppler factor for a sample of 105 sources, using data compiled 
from the literature. The data compiled by GPCM93 were used by Giiijosa & Daly (1996) to 
compute equipartition Doppler factors. These sets of inverse Compton Doppler factors and 
equipartition Doppler factors are used here to compute the Lorentz factors, viewing angles, and 
intrinsic properties of outflows from these sources. 

Updated estimates of both inverse Compton and equipartition Doppler factors for a sample 
of 100 sources are discussed in § |2.1| ; intrinsic brightness temperatures computed for both sets of 
Doppler factors are discussed in § |2.1.1| . The reliability of Doppler factors and the errors assigned 
to these estimates are discussed in j |2.2| and respectively. Estimates of the intrinsic luminosity 
density and intrinsic luminosity are examined in ^2.4| . The computations of the bulk Lorentz 
factor and the viewing angle are discussed in A subsample of sources with estimates of Doppler 



factors and apparent motion in the plane of the sky is introduced in §4.1, and a set of solutions for 



the bulk Lorentz factor and viewing angle are computed in § |4.2| for this subsample. In § [4.3| , these 
results are discussed in the context of orientation unified models. A general discussion follows in 
§1 



2. A SAMPLE OF DOPPLER FACTOR ESTIMATES 

2.1. Previous and Updated Doppler Factor Estimates 

GPCM93 assemble the relevant data to compute the inverse Compton Doppler factor for 
those AGN with VLBI core size data available in the literature circa 1986-92. The five BL Lacs 
without known redshifts are excluded from the sample examined here, since inverse Compton and 
equipartition Doppler factor estimates depend on redshift (see eqs. Al & A2). This leaves 100 
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sources in this sample which fall into the following classifications (GPCM93): 32 BL Lacerate 
objects (BL Lac), 53 core-dominated quasars (CDQ), 11 lobe-dominated quasars (LDQ), and 9 
radio galaxies (RG). The CDQs are further classified on the basis of their optical polarization 
into 24 core-dominated high-polarization quasars (CDHPQ), 22 core-dominated low-polarization 
quasars (CDLPQ), and 7 core-dominated quasars without polarization information (CDQ-NPI). 
A detailed discussion of these classifications can be found in GPCM93. 

Both the inverse Compton Doppler factor, 6ic, and the equipartition Doppler factor, 6eq, can 
be computed using the data compiled by GPCM93 (see the Appendix for the relevant formulae). 
By assuming that the observed radio frequencies and radio flux densities are the true peak radio 
frequencies and radio flux densities, GPCM93 use the data to compute 5ic- This approximation is 
made since the multi-frequency observations and analysis of the radio spectra necessary to obtain 
true peak frequencies and flux densities are available for only a small fraction of these sources (e.g. 
Marscher & Broderick 1985, Unwin et al. 1994). Giiijosa & Daly (1996) recompute 6ic taking into 
account corrections to the angular size and flux density neglected by GPCM93 (see Appendix), 
and compute 5eq with the same assumptions for this sample. 

Both 5eq and 6ic are recomputed in this study assuming a deceleration parameter of Qo = 0.05 
(i.e. Jlo = 0.1, r^A = 0) and h = 0.75, and using updated VLBI data for four sources (1101-1-384, 
0615-1-820, 1039-F811, & 1150-^812) from Xu et al. (1995); here is the mean mass density 
relative to the critical value, is the normalized cosmological constant, and Bubble's constant 
is parameterized in the usual way: Hq = 100 h km s~^ kpc~^. When computing 6eq and 6jc for 
each source, we assume an optically thin spectral index a = 0.75 {Si, cx v^") and a spherical 
geometry, as was done by GPCM93 and Giiijosa & Daly (1996) (see eqs. ^l] & |A2D . Note that 6eq 
is changed by only a few percent relative to the values computed assuming an Einstein de-Sitter 
universe (JIq = 1-0, f^A = 0) and h = 1, while 6ic is independent of these assumptions. Following 
GPCM93, upper bounds on angular sizes and lower bounds on redshifts are taken as detections in 
this section, but are treated as bounds when computing bulk Lorentz factors and viewing angles 



in §4.2. 



Figures |^a & |l]b show 5eq and 5ic respectively, computed as described in the Appendix, as 
functions of (l + z). In these figures and throughout, solid circles represent BL Lacs, solid diamonds 
represent CDHPQs, solid squares represent CDLPQs, solid triangles represent CDQ-NPIs, open 
diamonds represent LDQs, and open squares represent RGs. The error estimates for 5eq and 
6jc are discussed in § p.3| . The trend in redshift is easily understood in terms of the effect of 
Doppler boosting on the observed flux density and the flux limited nature of this sample. Only at 
low redshift are we able to observe those sources with lower Doppler factors, and the sample is 
dominated by higher Doppler factors at all redshifts. It is also apparent that most LDQs and RGs 
have Doppler factors less than 1, and most CDQs have Doppler factors greater than one. BL Lacs 
have a much wider range of Doppler factors, covering the entire range of the whole sample. 

Figure ^ shows the distribution of the equipartition Doppler factors, 6eq (dotted line), and 
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the inverse Compton Doppler factors, 5ic (dashed line), for all 100 sources in the sample; the 
horizontal axis is expressed in terms of log 5. Note that the distributions of 5eq and 5ic are similar, 
and both have a range of a few orders of magnitude. In addition, the peaks of both distributions 
are between 5 and 10. Figures |3|a-f show the distribution of 5eq (dotted line) and 5ic (dashed 
line) for each class of AGN described above. For BL Lacs, the distributions range over a few 
orders of magnitude in 5, with a peak around 5 ~ 5 (Figure ^a). CDHPQs have 5eq and 6ic that 
cover a smaller range of values than BL Lacs, and peak noticeably at about (5 ~ 10 (Figure |3|b), 
while the distribution of 5eq and 5ic for CDLPQs have a similar range as CDHPQs, and peak 
somewhere around 5 ~ 5 (Figure ^). CDQ-NPIs have the same range of 5eq and 5ic as CDHPQs 
and CDLPQs (Figure ^d). The distributions of 5eq and 5ic for LDQs are wider than those for 
RGs, but both classes peak at (5 < 1 (see Figures ^ & ^). 

The classification 3C216 as a LDQ should be examined closely since it has relatively large 
values for 5eq and 5ic (61 and 33 respectively). The separation of quasars into CDQ and LDQ was 
done by GPCM93 on the basis of the core dominance parameter computed using fluxes that were 
K-corrected to 5 Ghz. Recent observation by Reid et al. (1995) at 5 GHz give a core dominance 
parameter slightly greater than 1.0, and this would place 3C216 between the LDQs and CDQs in 
terms of the core dominance parameter. Most LDQs in this sample exhibit a triple morphology, 
consisting of a core with a flat radio spectrum and two lobes with steep radio spectra and have 
angular extents of about 10"-30" (e.g. Hooimeyer et al. 1992|) . 3C216 shows a triple morphology, 
but has an angular extent of only about 1.7". It could be argued that 3C216 should be reassigned 
as a CDQ, but in this paper we keep this source in the LDQ sample but make a special note of it. 



2.1.1. Brightness Temperatures 

Both the equipartition Doppler factor and inverse Compton Doppler factor, 5eq and 5ic, are 
used here to compute the intrinsic peak brightness temperature for each source in the sample 
(see the Appendix for relevant formulae). Histograms of the observed brightness temperature, 
intrinsic brightness temperature base on 5eq, and intrinsic brightness temperature base on 5ic 
{Tbo, TBi{eq), and TBi{IC)) are shown in Figure ^ for all 100 sources, and in Figures |5|a-f for each 
class of AGN. Table |l] lists the mean and the standard deviation of the mean of Teo, TBi{eq), and 
TBi{IC) for the full sample and each class. 

The observed peak brightness temperatures span a few orders of magnitude even within some 
individual classes of AGN. In contrast, there is a sharp peak in the distribution of the estimates 
of intrinsic brightness temperature, TBi{eq) and TBi{IC), and both are centered around about 
7.5 X 10^° K. These two estimates of the intrinsic brightness temperature agree well with each 
other, which is not surprising given the strong correlation between 5eq and 5ic described above 
and by Giiijosa & Daly (1996). 

It is interesting to note that the mean TBi{IC) for all classes (except CDQ-NPIs) are 
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consistent within la, while the LDQs and RGs have larger TBi{eq) than the other classes at 
about the 2a level. These two classes have values of R = 6eq/Sic less than one (Giiijosa &; Daly 
1996), which could possibly be due to a systematic underestimate of 6eg by about 75%. Such an 
underestimate of 5eq could account for the larger Tsiieq) found for LDQs and RGs. 

The standard deviations of Tsiieq) and TBi{IC) are 0.27 x lO^^ K and 0.33 x 10^^ K 
respectively. It should be noted that TBi{eq) oc S^piy^pO'^'^ and TsiilC) oc S^pU~p-'^ 6'^^''^ which 
are weak dependences on observable, much weaker than the dependences of Tbo on the same 
parameters. If the range of values for the observables that go into Tsiieq) and Tbi^IC) are not 
too large, then their narrow distributions can be explained by weak dependences on observable 
parameters. The range of 9d contributes the most to the range of TBi{eq) and TBi{IC), since Uop 
is fixed at discrete values corresponding to the frequencies of observations, and Sgp is virtually 
insignificant in computing TBi{eq) and Tbi^IC). If the narrow ranges of TBi{eq) and Tbi^IC) are 
due to some physical cause and not the range in 9d, then randomly reassigning the angular sizes 
used to compute TBi{eq) and TBi{IC) should widen the distributions of each. When the angles are 
reassigned randomly, the means are quite similar (0.85 x 10^^ K and 0.79 x 10^^ K respectively), 
and the standard deviations are 0.26 x 10^^ K and 0.43 x 10^^ K respectively. Each width is 
not significantly affected by randomly reassigning angular sizes, which suggests that the narrow 
range of TBi{eq) and TBi{IC) obtained here is primarily due to the range of observables, and are 
probably not due to the physical state of the sources. 

Readhead (1994) points out that if the intrinsic brightness temperatures were cutoff due to 
the inverse Compton catastrophe (i.e. the radiation lifetime plays a role in the observed cutoff), 
one would expect the intrinsic brightness temperatures of compact radio sources to cluster near 
]^q11.5 Both estimates of intrinsic brightness temperatures for this sample peak below 10^^'^ K 
which agrees with the results of Readhead (1994). This suggests that the inverse Compton 
catastrophe does not play a role in the observed brightness temperature distribution, or if the 
peaks are set by the inverse Compton catastrophe, a factor of about 4 has not been accounted for 
in computing the unbeamed brightness temperatures. 

2.2. Reliability of 5eq and 5ic as Estimators of the Doppler Factor 

Previous results indicate that 5eq — 5/c for this sample (Giiijosa & Daly 1996). To further 
test this conclusion, the effect of randomly reassigning one of the observable quantities used to 
compute 5eq and 5ic is examined. Since both Doppler factors depend strongly on angular size 
{5eq oc 0^^'^ and Sjc oc 6*^^'^), the angular sizes used to compute 6eq and 5ic were randomly 
reassigned between sources within each class, and the mean values are shown in Table ^. Note 
that the means of 5eq and 5ic tend to systematically increase from their true values for all classes 
when the angular sizes are reassigned. 

The fractional uncertainties in the mean values of 5eq and 5ic are of order 20% to 30%, so 
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that the mean 5eq and 5ic are significant at about the 3a to 5a level (Giiijosa & Daly 1996). 
In contrast, the "random" mean values are significant at 2a or less (i.e. their uncertainties are 
greater than 50%, Table |2|). This suggests that these Doppler factor estimate are reliable. 



2.3. Errors for the Doppler Factor Estimates 

The main source of error for both Doppler factor estimates is the assumption that the radio 
flux density and frequency used to compute the Doppler factor are the peak radio flux density and 
peak frequency. Thus, an estimate of errors for this data set should focus on systematic errors 
instead of statistical ones. A dependence of observed angular sizes on frequency occurs if the 
synchrotron-emitting plasma has gradients in magnetic field and particle density (Marscher 1977, 
1987), but since the model considered here is that of a uniform sphere, the errors associated with 
this frequency dependence of angular size are neglected. 

The computed errors are found by assuming all sources have the same fractional errors on 
the peak frequency, Uop] this uncertainty is likely of order 100% and dominates all other sources 
of error, as discussed in Giiijosa &: Daly (1996). It is assumed that the flux density has an error 
which is related to the spectral index between the observed frequency and the peak frequency, 
Oop, where Sop = So{i'op/^o)~"'°^' and So and Uo are the observed flux density and frequency; Uop 
should not to be confused with the optically thin spectral index, a. In terms of VLBI observed 
properties, 5eq oc Sl^i'^'^9'^^'^ and 6ic oc Sl^u'^'^O'^^'^, assuming a = 0.75. These two Doppler 
factor estimates can be related to the ratio of the peak frequency to the observed frequency, 
6eq OC {i'op/j^o)~^'^°'°''~'^'^ and 6ic oc (z^op/'^o)"^'*^"""""^'^- One finds that 

^ = il.laop + 2.3)(^] and ^ = (1.0a„p + 1.3) f (1) 

deq \ fop J Oic \ Vop J 

where aSeq-, a Sic-, and auop are the error estimates for (5eg, and Uop- 

An estimate of avop/^op is needed to assign errors for the Doppler factor estimates using 
equation (|^). A reasonable value of avop/^op is that which gives a reduced of 1.0 when fitting 
R = 5eq/5ic to unity. 

Since R oc Sl:^v-^-^e^^-'^ for a = 0.75, it can be shown that R oc (z^op/j^o)"°-^"°''"^-°. The 
error on the ratio R, denoted aR, can be related to avop/^op by 

Note that if auop/fop and Oop are the same for all sources, then aR/R is the same for all sources. 

The value of aR/R = 0.61 gives a reduced of 1.0 when fitting i? to a constant equal to 
unity for all sources. A value of Oop = 0.4 is assumed for all sources since radio sources with 
observed optically thin spectral indices for their cores have spectral indices in the range of about 
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0.5 to 1, and the spectral index must go to zero towards the pea-k. A lower value of cxop 

could be 

used, and this will decrease the errors slightly on 5eg and 5ic- For Oop = 0.4 and aR/R = 0.61, 
equation (||) gives (JVop/i^op = 0.59, and one finds that a5(,q/Seq = 1-6 and odic/^ic = 1-0 using 
equation (|l|). These error estimates for 5eq and 5ic are adopted throughout this paper. 

Some objects in this sample have been observed at multiple frequencies, and spectra are 
available in the literature. The observed frequencies used to compute the Doppler factors for 
eight of the sources studied here are compared to the peak frequencies of the integrated spectra 
available in Xu et al. (1995), in order to check if our estimate of (Jfop/^op is reasonable. It should 
be noted that the integrated spectra include components besides the cores, and that the core 
component is used to determine the Doppler factor (see §^). Here we approximate that the peak 
of the integrated spectra corresponds to the peak of the core. On this basis, four sources have 
ouoplvfyp between and 0.5, and the other four have ovopjvop between 0.5 and 4, which would 
indicate that the value of avoplvop = 0.59 adopted for the whole sample is quite reasonable. 



2.4. Luminosity Densities and Luminosities 

The uncorrected luminosity density is given by 

=Maor)2(l + ^)5o, (3) 

where So is the observed flux density, z is the source redshift, (ao^) is the coordinate distance to 
the source, and v^n is frequency as would be measured in the local Hubble frame of the source 
without correcting for Doppler boosting (f„„ = z^o(l + z) where Vq is the observed frequency). 
The flux density corrected for Doppler boosting of a uniform spherical source, Sj, is given by 
Si = So/ 6^, where 5 is the Doppler factor. The corrected flux density is the flux density that 
would be received by an observer moving with the same velocity with respect to the Hubble flow 
as the radiating source. Substituting Si for So in equation (^) gives the Doppler corrected, or 
intrinsic, luminosity density: 

1 + z 

Lyi = A-K{aorf—^So (4) 

where Ui = Uo{l + z)/5. Substituting both estimates of the Doppler factor d^q and 5ic into 
equation (Q) gives two sets of intrinsic luminosity densities which shall be referred to as Lyi{eq) 
and Lyi{IC) respectively. 

It should be noted that equations @ and (^ are the luminosity densities at frequencies that 
differ from the observing frequencies, which are = ^o{^ + and vi = Vo{^ + z)/^ respectively. 
Thus, the uncorrected luminosity can be approximated as L^n ~ L^^^Uoil + z), and the intrinsic 
luminosity can be approximated as 



Li L^ii^oil + z)/6. 



(5) 
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Note that luminosities computed in this manner are crude approximations that may be 
systematically biased by the approximated peak frequencies, yet still give order of magnitude 
estimates. 

Figures § & |^ show Lj^^^ and Lun as functions of (1 + z). The flux- limited nature of this 
sample is evident in these two figures, and it is quite obvious that selection effects play a role in 
the composition of this sample. The intrinsic luminosity densities, L^i{eq), are plotted versus the 
rest frame frequency i'i{eq) in Figure ^ and Li{eq) versus i^i{eq) is plotted in Figure By fitting 
all the data plotted in Figures § and one finds that L^i{eq) oc fj(eg)^ '^^'^'^ with a reduced 
equal to 4.6 and a correlation coefficient r = 0.91, and Li{eq) oc i/j(eg)^'^^'^'^ with a reduced 
equal to 4.0 and a correlation coefficient r = 0.96. Both of these power law fits are very significant 
since r is close to unity. Figures IC and ^ show Lyi{eq) and Li{eq) respectively as functions of 



(1 + z); fitting all these data, one finds that Lyi{eq) oc (1 + z)"^'^^^'^ with a reduced equal to 

1.3 and a correlation coefficient r = 0.16, and Li{eq) oc (1 + z)^-^^"^-"^ with a reduced equal to 

1.4 and a correlation coefficient r = 0.04. Similar results are found for Lui(IC) and Li{IC). 

Figures § & ^ show a trend in Lyi{eq) and Li{eq) with Vi{eq) that can be explained by the 
strong dependence of these quantities on the Doppler factor and a systematic effect introduced 



by assuming the observed frequency is the peak frequency. Equations (0) & (A2) indicate that 



oc 5eq OC v^p and that Viieq) oc i^op^eq oc v^p , for a — 0.75. Since it is assumed that 
i^op = i^o when computing 6eq, the offset of Vo/^op from unity will cause Lyi{eq) to be offset by a 
factor of (fo/fop)^'^ and I'iieq) to be offset by a factor of (vo/^op)^'^ from their actual values. One 
would expect to find a systematic effect that would increase the range of L^i{eq) and i^i{eq) such 
that Lyi{eq) oc i'i{eq)'^'^ . In fact, one obtains Lyi{eq) oc Vi{eq)'^'^^^'^ by fitting the data in Figure]^. 
A similar calculation gives the prediction that Li(eq) oc fj(eg)'^'^, and the fits of the data in Figure 
^ give Li{eq) oc Vi{eq)^'^^^'^ . The spread over such a large range of values for Lyi{eq), Li{eq), and 
i'i{eq) can thus be explained by this systematic effect. On the other hand, the errors in Figures ^ 
& ^ are quite large, and all the Lyi{eq) and Li{eq) are still consistent with each other. 

It is interesting to note that Iu{eq) oc Vi{eqY''^^^'^ and 1^(10) oc fj(/C)^'*^^°'^ for this sample, 
where I,^ oc Lyi/O^'- This dependence would be expected if these sources were intrinsically 
similar, have black a body spectrum, and are observed in the Rayleigh-Jeans size of the black 
body spectrum. It should be noted that the small dispersion in intrinsic brightness temperatures 
described above (§ |2.l| ) and the fact that I^, is approximately proportional to v'^ go hand in hand. 



Table ^ lists the mean and standard deviation of the mean for the logarithmic values of 
Lyo, Lo, Lt^i{eq), Li{eq), Lyi{IC), and Li{IC), for each class of AGN. Comparing Lyi{eq) to 
Lyi{IC) and Li{eq) to Li(IC), agreement for each class of source is found within about 2a. 
Ranges of intrinsic luminosity densities and intrinsic luminosities for all classes of AGN are 
^032.4to33.ig^gs g-iHz-i and W^'^ ^^ergs s'^. 
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3. DISENTANGLING RELATIVISTIC EFFECTS IN RADIO-LOUD AGN 

The effects of relativistic motion can be described by two quantities: the bulk Lorentz factor 
of the outflow, r, and the angle subtended by the outflow direction and the line of sight to the 
observer, (j). Note that the bulk Lorentz factor is determined by the speed of the outflow relative 



to the speed of light, /3, by the familiar equation T = 1/v 1 — /?^- 

Two observable quantities which are combinations of T and (j) are the Doppler factor, 5, and 
the apparent speed projected onto the plane of the sky relative to the speed of light, (3app- These 
are expressed in terms of (3 and (j) as 



(6) 



1 — /? cos (j) ' 
and 

/?sin0 
^"PP 1-/3COS,/.- 

Equations (|6|) and (^) can be used to solve for F and in terms of 5 and (iapp-, and it can be shown 
that 

+ + l 

r = (8) 

and 

yapp I ^ ^ 

(see GPCM93; Daly, Guerra, & Giiijosa 1996). Thus, F and (/> can be estimated for relativistic 
outflows in AGN, by combining observational estimates of 5 and Papp- Here 5eq and Sjc are used 
to estimate 5. 



3.1. Apparent Bulk Flow Versus Pattern Flow 

It has been suggested that the pattern velocities, the velocities of radio features observed in 
VLBI proper motion studies, are greater than the bulk velocities, the velocities responsible for the 
observed Doppler boosting (see, for example, Lind & Blandford 1985). This issue was addressed 
by GPCM93, who conclude that the data are consistent with pattern velocities equaling bulk 
velocities for the sources in their sample with observed superluminal motion. 

Vermeulen & Cohen (1994, hereafter VC94) examine the equivalence of pattern velocities and 
bulk velocities in detail using Monte Carlo simulations of a simple model where the angle between 
the direction of outflow and the line of sight is allowed to vary randomly, but the bulk Lorentz 
factor is the same in every jet. They find that the pattern Lorentz factor, Fp, would have to be 
about twice the bulk Lorentz factor, Ff,, in order to be consistent with the data from the GPCM93 
sample. They point out, however, that for a broad distribution of bulk Lorentz factors the data 
are easily consistent with the Fp ~ F;,. 
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It is assumed, here and throughout, that Tp ~ r^, and any departures from this equahty are 
of order unity. 



4. ESTIMATES OF BULK LORENTZ FACTORS AND VIEWING ANGLES 

4.1. A Subsample with Papp Compiled 

VC94 compile a hst of 66 radio sources with multi-epoch VLBI observations of the internal 
proper motions. Table 1 in VC94 lists Papp for the features in these radio sources. There are 43 
sources that overlap the GPCM93 and the VC94 samples, and these sources are examined below 
in order to estimate T and (j) (§^]^); preliminary results were presented by Daly, Guerra, & Giiijosa 
(1996). 

Table § lists the 43 sources that overlap the GPCM93 and VC94 samples. Column (1) lists 
the lAU Name, column (2) lists the common name if it exists. Column (3) gives the redshift listed 
by VC94, Column (4) lists the Papp estimates from VC94, except for 0108+388 where the estimate 
is based on more recent observations by Taylor, Readhead, & Pearson (1996). When multiple 
components are listed by VC94, the weighted mean is taken for Papp- Figure |l2| shows Papp versus 
(1 + z) for the 43 sources in the overlapping sample. Columns (5) and (6) give estimates of 6eq 



and 6iCi respectively, which are discussed in §2.1 above. 



4.2. Values for F and cp 

The values of Papp, Seq, and 5ic listed in Table ^ can be used with equations (^) and (^) 
to produce two sets of estimates for T and (p. In the cases where only bounds on Papp or 5 are 
available, bounds on F and (p are computed. 

Table ^ lists the solutions for F and (p using the equipartition Doppler factor, called Feg and 
(peq, and using the inverse Compton Doppler factor, called F/c and (pic- Columns (1) and (2) 
list the lAU name and the common name respectively, Column (3) gives the redshift listed by 
VC94, Columns (4) and (5) list Fgg and (p^g respectively, and columns (6) and (7) list Tjc and (pjc 
respectively. Both sets of solutions for F and (p agree within errors for any given source in Table |. 

Table ^ lists the median values of T^q, (peq, Fjc, and (pjc for each class of source, excluding 
those values which are upper or lower bounds. The median values of the equipartition set {Teq and 
(peq) and the inverse Compton set {Tjc and (pic) of estimates agree well, except for the median 
Teq and Tjc of the LDQs which are a factor of about two greater in the equipartition case. Note 
that excluding 3C216 from the LDQs would not significantly change the median values for this 
category. 

Figures |l^ a & 13b show Teq versus (peq with and without errors respectively, and Figures |l^ a 
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Sz |14| b show Tjc versus (j)jc with and without errors respectively. The figures without errors are 
shown so that the symbols are more visible to the reader and bounds on F and cp are denoted by 
arrows in these figures. Figures p^ a & [T5| b show (j)f,g and Fg^ as functions of (1 + z) respectively, 
and Figures 16a & 16b show (pjc and Tjc as functions of (1 + z) respectively. 

The median values for both sets of estimates of F and (p make it possible to compare different 
classes of AGN in the context of orientation unified models (see discussion in §[4.3|). A detailed 
discussion of each set of F and (j) estimates follows in §§4.2. l| & [4.2.2. 



4-2.1. Solutions Using the Equipartition Doppler Factor 

The eight BL Lacs studied here span a relatively large range of Fe^; one source has a value of 
Tgq = 1.4 lb 0.2, and two sources have Fgg between 40 and 100, although these two sources have 
large errors on Fgg. The remaining five BL Lacs have Teq between 3 and 6. The values of (peg for 
BL Lacs span a range from about 0° to about 50°, with six of the eight BL Lacs having (peg ^ 20°. 

Seven of the eight CDHPQs have (peg ^ 12°, and at least five of these sources have (peg ^ 6°. 
Three out of eight CDHPQs have T^g between 3 and 10, and two CDHPQs have T^q between 20 



and 40. Three sources have limits on Fgg and (peq instead of values (see §2J.). Two of these sources 
with bounds, 1156+295 (4C 29.45) and 2230+114 (CTA 102), stand out in Figures |l|a,b. The 
bounds on Teq for these two sources, Teq < 405 for 1156+295 and Teq > 2.8 for 2230+114, still 
allow these sources to have Lorentz factor similar to other CDHPQs. 

Five CDLPQs have Teq between 5 and 10, and three have Teq between 15 and 30. The other 
two sources only have bounds on Teq. Eight out of ten CDLPQs have (peq ^ 14°. 

The three CDQ-NPIs all have Teq < 9, and have quite different values for (peq. The sources in 
this small category in reality belong to the CDHPQs or CDLPQs, and may contain both types. 

Five LDQs have Teq between 10 and 50. These values of Teq are larger than the typical 
values for the other classes of AGN studied here. One source has Teq less than 10 (4C 21.35, 
Teq = 2.8 ± 2.9), while another has a lower bound Teq > 3. Only one LDQ has (peq less than 15° 
(3C216), while the other six LDQs have (peg between 15° and 41°. LDQs have the larger values for 
(peq than most of the BL Lacs, CDQs (HP, LP, and NPI) studied here. 

RGs have the largest values and widest range of (peg of all the categories. Four out of seven 
RGs have 45° < (peq < 135°, which suggests that these sources have outflows that typically lie 
close to the plane of the sky (see §|4.3|). RGs also have the lowest Teq of all the categories, with 
at least four RGs having Teq < 2. Note that the radio galaxies in this sample are compact steep 
spectrum sources or FRI; there are no FRII radio galaxies in this sample. 
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4-2.2. Solutions Using the Inverse Compton Doppler Factor 

The results for Tjc and (pjc are quite similar to those for Teq and (j)eq, and agree well within 
errors. This agreement is not surprising considering the correlation and agreement found between 



6ic and 6eq (see §2.1). 



Five out of eight BL Lacs have Tic between 3 and 5, while the total range for all BL Lacs is 
between 1.5 and 60. Seven BL Lacs have (pjc < 20°. Five out of eight CDHPQs have Tjc between 
8 and 16, while the other three have various upper and lower bounds. These same five CDHPQs 
all have (pic < 10°. Six out of ten CDLPQs have Tic between 3 and 11, and two have Tic around 
20. Eight CDLPQs have (pic < 14°. Five out of seven LDQs have Tic between 10 and 30, and 
five LDQs have (pic between 15° and 45°. As in the equipartition case only one LDQ has (pic less 
than 15°. Three out of seven RGs have (pic between 45° and 90°, and three RGs have (pic > 90°. 
At least six RGs have Peg < 2. 



4.3. Implications for Orientation Unified Models of AGN 



The results presented §4.2 are consistent with an orientation unification scheme for radio-loud 
AGN (e.g. Antonucci 1993). Common orientation models have a torus or disk of material that 
absorbs or obscures the broad absorption lines and other non-stellar radiation emitted from the 
nuclear region of the AGN, and the classification as a radio galaxy occurs if the torus intersects 
the line of sight from the nuclear region to the observer. The radio jet axis in radio-loud AGN is 
thought to be more or less perpendicular to the plane of the torus or disk. In such models (e.g. 
Barthel 1989), RGs, LDQs, and CDQs, are viewed with different angles with respect to the jet 
axis {(p, which is referred to as the viewing angle): radio galaxies have the largest angles {(p ^ 45°), 
LDQs have smaller angles (20° ^ (p ^ 45°), and CDQs have the smallest angles {(p ^ 20°). 

The RGs, LDQs, and CDQs in overlap of the GPCM93 and VC94 samples have estimates of 
(p consistent with this orientation unification scheme. The RGs in this sample have a median (p^q 
of 110° and a median (pic of 81°, which suggests that they have jets that lie close to the plane 
of the sky. LDQs typically have a median (peg of 26° and a median (pic of 25°. CDLPQs have a 
median (peg of 7° and a median (pic of 6°, while CDHPQs have a median (peg of 4° and a median 
(pic of 3°. The median values for CDQ-NPIs are insignificant since only two sources are used to 
compute them. It is interesting to note that CDHPQs and CDLPQs have similar median (peg and 
(pic- This suggests that it is intrinsic differences between these two classes of AGN (not viewing 
angles) that determine the amount of optical polarization observed. 

The BL Lacs in this sample cover a wide range of F and (p, but have viewing angles typically 
less than 45°. An extension of the orientation model described above assigns BL Lacs to the 



parent population of FRI (edge-darkened, Fanaroff & Riley 1974) RGs, where BL Lacs are seen 



at smaller viewing angles and radio selected BL Lacs have smaller angles than X-ray selected BL 
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Lacs (see Padovani &: Urry 1992). The sample of BL Lacs examined here are almost all radio 
selected. Mkn 421 (1101+384) is the only X-ray selected BL Lac in this sample, but does not 
stand out against the other BL Lacs in terms of its T or 4> estimates. 

The BL Lacs in this sample have a median 0eg of 12° and a median (pjc of 14°. These 
results are consistent with the those found by Kollgaard et al. (1996), who estimate the viewing 
angles for samples of radio selected BL Lacs, X-ray selected BL Lacs, and RGs by comparing the 
core enhancement relative to the more diffuse radio emission for each class of AGN. Assuming a 
constant bulk Lorentz factor, they find that their samples of radio selected BL Lacs, X-ray selected 
BL Lacs, and RGs have average viewing angles of 10°, 20°, and 60° respectively. 

CDQs in this sample are centered around a L ~ 9. If it is assumed that the distribution of 
bulk Lorentz factors for all classes of radio- loud AGN is centered around F ~ 9, then one might 
naively expect the mean or median F for each of class of AGN to similar. Selection biases likely 
play a role in the composition of the samples studied, and systematic effects can affect the values 
that can be obtained from the data. 

LDQs have median T^q and Tjc that are larger than those for the CDQs, and these two 
estimates are different almost by a factor of two, yet the mean ratio of i? = Seq/Sjc is only about 
0.76 for LDQs. This could be caused by the sensitivity of the estimates of F to the precise values of 
Papp and 6 when these to parameters are close to one (see eq. which is the case for this sample 
of LDQs. A decrease of 6 away from one will cause the corresponding F estimate to increase. 
If the 6 estimates of LDQs in this sample are systematically underestimated, then a noticeable 
systematic increase in the F will result, which could be caused by an overestimate of the peak 



frequency (see eqs. A2 and Al). It is not clear whether this is one of the causes of the larger F for 
LDQs. 

The lower F for RGs compared other classes of AGN can be explained as follows. RGs will 
appear in this sample if there cores are bright enough to have VLBI core size data, which occurs if 
the sources are nearby (like M87 and NGC 6251) or have very bright compact cores. In particular, 
FRII radio galaxies are missing from this sample since a majority of their cores are not significantly 
Doppler boosted and they make difficult targets for VLBI. In the unification scheme described 
above, RGs have viewing angles around 90°. Approximating (p ~ 90°, one finds that 6 ~ F~^ 
from equation (P). Thus, F ~ and those RGs with smaller F will appear in this sample since 
these RGs will have larger 6 and brighter cores. Even if the true distribution of F for RGs peaked 
at around 10, we may only be able to observe RGs with F from 1 to 2. It should be noted that 
RGs in this sample have values of F near unity which are in agreement with other evidence that 
Doppler boosting is not significant in these sources (e.g. Taylor, Readhead, & Pearson 1996). 

There are two sources that seem to have viewing angles greater than 90° based on equipartition 
and inverse Compton Doppler factors. One, 2352+495, has a lower bound on the viewing angle, 
the other, 0710+439, is greater than 90° at around 3a. This would suggest that some RGs have 
observed outflows that are moving away from the observer, and that if there is an outflow moving 
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toward the observer it must be intrinsically fainter. 

BL Lacs have lower values of T than the CDQs, with a median Teg of 4.2 and a median Tjc 
of 5. Since the BL Lacs studied in this paper are almost exclusively radio selected, it would be 
interesting to compare T and (p estimates for a sample of X-ray selected BL Lacs. Kollgaard et al. 
(1996) find that the core enhancements they observe for BL Lacs and RGs, require that T > 4.5, 
assuming a constant bulk Lorentz factor for all these sources. The estimates of F found here are 
barely consistent with T > 4.5, but a distribution of Lorentz factors may loosen the constraints 
placed by Kollgaard et al. (1996). 



5. DISCUSSION 



The two estimates of the Doppler factor for the sample of radio-loud AGN discussed in § p.l 
agree with each other on average and provide a means of estimating the intrinsic properties of 
these sources (e.g. brightness temperatures, luminosities) and the parameters that describe the 
kinematics of the relativistic bulk flow of the radio emitting plasma. Although the results here 
are based on rough estimates of the Doppler factor, the average or typical properties of different 
classes of AGN can be compared. 

Care should be taken when using the estimates of Doppler factors and intrinsic properties 
computed here since systematic errors contribute most of the uncertainty. These systematic errors 
arise from the simplification that the observed frequencies correspond to the peak of the core 
emission. Another concern with this sample is completeness. The sample examined here was 
defined by GPCM93 as those which had VLBI core size data in the literature at that point in 
time. The sources that enter this sample will have VLBI cores with higher surface brightness than 
those excluded from the sample; this will favor the brighter (more highly Doppler boosted) and 
more compact (smaller angular size) cores. GPCM93 warn that this sample is not complete in any 
statistical sense, and that biases may have a large effect on the content of this sample. 

Sources in the first and second Caltech-Jodrell Bank VLBI surveys (CJl and CJ2) form 
a complete flux-limited sample of almost 400 compact radio sources with VLBI measured core 
sizes (see Henstock et al. 1995, Xu et al. 1995) and some multi-frequency data. Doppler factor 
estimates for sources from these samples are currently under investigation, including the use of 
spectra to estimate the peak frequency of core components in order to more accurately compute 



Doppler factors ( Guerra fc Daly 1997 ). Large samples such as these could be used to estimate the 
Doppler factors and intrinsic properties of different classes of AGN with greater confidence and 
statistical weight. 

The intrinsic brightness temperatures for this sample are estimated using two sets of Doppler 
factor estimates, 6eq and 5ic, and both histograms of these values are centered around 7.5 x lO^'^ K 
(see § p.lD . The intrinsic brightness temperatures are similar for different classes of AGN, despite 
the fact that the observed brightness temperatures are different for different classes of AGN. The 
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center of the distribution of intrinsic brightness temperatures is a factor of 4 lower than what 
would be expected if the "inverse Compton catastrophe" played a major role in the observed 
brightness temperature cutoff, as noted previously by Readhead (1994). This suggests that either 
a factor of about 4 has not been accounted for in these estimates (e.g. a geometric factor), or 
that a physical process other than the "inverse Compton catastrophe" is limiting the intrinsic 
brightness temperature. 

Estimates of the intrinsic luminosity density and luminosity depend strongly on the Doppler 
factor {Li,i oc and Li oc (5~^), and systematic errors are quite large for these estimates. The 
large scatter of Li^i and Lj is a cause for concern since it can be attributed to systematic errors 
introduced by assuming the observed frequencies correspond to the actual peak frequencies of these 
sources. In any event, these estimates can be compared to the Eddington luminosity for a supposed 
massive central object. The Eddington luminosity is Le = 1.3 x 10*^^ Mg ergs s~^ where Mq is 
the mass of the central object in units of 10^ Mq which gives a range from 10^"^ to 10^6 ergs s 
for central objects with Mq = 1 to 100. The mean intrinsic luminosity for the sources examined 
here is around 10'*^ to 10"^^ ergs s~^, while only five objects (3 BL Lacs, 1 LDQ, and 1 RG) have 
intrinsic luminosities greater than 10^^ ergs s~^ by more than la. This is consistent with the radio 
luminosity produced being 0.1 to 1 percent of the Eddington luminosity for a central compact 
object with mass 10*^ to 10^ Mq. 

A subsample of 43 sources have observed proper motions in the plane of the sky and were 
used by VC94 and others to compute Papp- The values of Papp for each of these sources can 
be combined with the Doppler factor estimates to solve for the bulk Lorentz factor, T, and the 
viewing angle, (j) (see §|3|). The estimates using the data compiled here are consistent with the 
orientations unification scheme for AGN (discussed in § |4.3| ), where CDQs have the smallest (j) 
(< 15°), LDQs have larger (p (15° < 4> < 45°), and RGs have the largest 4> (> 45°); note that the 
RGs in this sample are compact steep spectrum sources or FRI. 

Selection biases and systematic errors may play a role in these estimates of T and 4>, in 
particular for LDQs and RGs (see §^|^). The RGs in this sample have lower T than the other 
classes of AGN, which is consistent with the parent population having outflows that lie close to 
the plane of the sky (i.e. (p ~ 90°) and the flux- limited nature of the sample studied. The LDQs 
are in a regime where small systematic errors in 6 can strongly affect T estimates. 

The estimates of (j) and T for the BL Lacs in this sample are consistent with the results of 



Kollgaard et al. (1996) for radio selected BL Lacs (see §0). Their estimates of (j) are based on 



the core enhancement of radio selected BL Lacs, X-ray selected BL Lacs, and RGs. A sample of 



X-ray selected BL Lacs with estimates of (j) and F computed as in §4.2 would give a consistency 
check for both methods. 

Marscher (1977) discusses the effect of gradients in magnetic field and density with radius 
on the relationships between the observed quantities and the important physical parameters 
in a nonuniform source. The effect on 6eq and 5ic of including inhomogeneities is constant 
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multiplicative factor from 0.7 to 1.4, and the the functional dependence on observable quantities is 
identical to the homogeneous case; the form of the inhomogeneity and the parameter choices are 
discussed by Marscher (1977). It should be noted that the intrinsic brightness temperatures will 
be changed by the same constant factor, and that inhomogeneities do not account for the factor 
of 4 in brightness temperature that would allow the inverse Compton catastrophe to play a major 
role in the brightness temperature distribution ( |2.1| ). 



The spherical model for the radio cores of AGN, which is used throughout this paper to 
compute Doppler factors, may not be as realistic as those models which consider a jet-like 
geometry. The expression for 5jc in the jet-like case, which is derived by GPCM93, is 

Sicijet) = [6ic{sph)f-^^''^/^^-^'''^ , (10) 



where 6ic{sph) is the expression given by equation (Al). For a value of a = 0.75, 
^icU^t) = 5ic{sphY '^ which is a power law only slightly greater than unity. The expression for 
5eq in the jet-like case is 

/ n \ l/(9+2a) 

6eq{jet) = [5e,(.p/.)](i3+2")/(9+2-) (sin,^)V(9+2.) (^1^ ^ (11) 



where 5eq{sph) is given by equation (A2), cj) is the angle subtended by the direction of the outflow 



and the line of sight, and {9a/9b) is the ratio of the angular sizes of the major and minor axes 
(see Appendix for the derivation of this equation). For a = 0.75, 6eq{jet) oc (sin^)''-^^, and 
^eqij^t) oc {Oa/Ob)^'^'^ , which are very weak dependences. The relation between the jet-like and 
spherical cases of the equipartition Doppler factor is 5eq{jet) tx 5eq{sphY''^ which is a power-law 
close to unity. The expressions for inverse Compton and equipartition Doppler factors in the 
jet-like case appear to be quite similar to those in the spherical case. Thus, approximating a 
spherical geometry for simplicity can still provide useful Doppler factor estimates even if the true 
geometry of these outflows is more jet-like. 
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Relevant Formulae 



For the simple case of a uniform sphere where the synchrotron-emitting particles have a 
power-law energy distribution and move through a tangled homogeneous magnetic field, the 
Doppler factor needed to reconcile the predicted and observed X-ray fluxes (Marscher 1987) is 

l/(4+2a) 



6ic{sph) = f{a)S„ 



In {ub/i^, 



op J 



SxGa 



Vo-p 



(1+^) 



(Al) 
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where z is the redshift, Sx is the observed X-ray flux density (in Jy) at frequency (keV), Qd is 
the angular diameter of the sources (in mas), z/^ is the synchrotron high-frequency cutoff (assumed 
to be 10^ GHz), a is the optically thin spectral index (where the radio flux density Sy oc v~^^^ i^op 
is the observed frequency of the radio peak (in GHz), and /(a) = 0.08a + 0.14 . The radio flux 



density in equation (Al), Sm (Jy), is the value obtained by extrapolating the optically thin flux 
density to the observed peak at fop- In equation (Al) and throughout this work, 9ii = l.SOpwHM 



and 0FWHM = V^adb from VLBI observations, where 9a and 9b are the angular sizes of the major 
and minor axes respectively (see Giiijosa & Daly 1996). 

The equipartition Doppler factor is estimated by assuming the energy density in magnetic 
field equals the energy density in synchrotron-emitting relativistic particles (Readhead 1994). In 
the notation used here, 

/r n34r i2 \ l/(13+2o) 

6e,{sph)=[[w^F{a)] [Ah/y{z)] {l + z)''>+'"SlX"'^op-'') , (A2) 

where Sop is the observed peak flux density and Uop is the observed frequency of the peak. 
The solutions for F{a) are given by Scott & Readhead (1977), and the solution of interest 
here is F(0.75) = 3.4 since a = 0.75 is assumed throughout. The function y{z) = Hoaor{z)/c 
(Peebles 1993) is a dimensionless function of Qo, ^A, and z, and contains the dependence on the 
coordinate distance to the source. In equation ( [A2[ ), the Hubble's constant is parameterized as 
Hq = 100 h km s^^ kpc~^^. The dependence of 6eq on cosmology is very weak. For a = 0.75, 
6eq OC {aor{z))~^'^'^ oc {h/y{z))^'^'^ Throughout this paper, Qo = 0.1, Qa = 0, and h = 0.75 are 
assumed, where the functional form of y{z) is that derived from equation (13.36) of Peebles (1993). 
Note that the results of Giiijosa & Daly (1996) were obtained assuming an Einstein-de Sitter 
cosmology (Oo = 1-0, = 0) with h = 1. 



The expression for 6eq{jet) (eq. |llD can be derived by equating equipartition magnetic field, 
Beq, to the synchrotron self-absorbed magnetic field, BssA, as is done when deriving 6eq{sph) 
(Readhead 1994). The equipartition magnetic field is given by 

B,q = (2 C„e h{a) [1 - exp(-rp)] Ip i// r'f/' (A3) 

where Cme and ii are defined by Leahy (1991), Tp is the optical depth of the peak of the spectrum 
(see Scott & Readhead 1977), is the frequency of the peak as measured in the frame of the 
synchrotron-emitting plasma, Ip is the peak specific intensity as measured in the plasma frame, I 
is the line of sight depth as measured in the plasma frame, and a is the optically thin spectral 
index (/jy oc z/~"). The synchrotron self-absorbed magnetic field is given by 



BssA 



"^^^"^,Ml-exp(-r,)]-i)%,-^.,^ (A4) 



.6 C2(a) 

where ci(a) and 02(0) are tabulated by Marscher (1977). To simplify these expression we define 

Geq{a) = (2 Cme h{a) Tp^ [1 " exp(-rp)] )2/7 (A5) 
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so that Beq = Geq{a){Ip Vp" I '^)'^/'^ , and 

,6 C2(a) 



eq — eqy-^jy^p "^p 

" 2 n ^^i' ^ M-i 



G55A(a) = ( J ^ r2 [1 - exp(-Tp)]-i)' (A6) 



'2 ,,5 



SO that BssA = GssA{a) Ip 

The equation for 6eq{sph) is derived by setting B^q = BssA, Ip = Sop (1 + 2:)'^ S^^ 6^"^, 
^p = ^op (1 + -z) '5""'^, ^ = Oo^'C-z) (1 + z)~^, and solving for 5. It should be noted that a spherical 
geometry is assumed when using these expressions for Ip, fp, and I (see GPCM93). The result, 

5,q{sphf^+^^ = Geq{af GssA{a)-' [aor{z)]'^ (1 + z)''^^'' S^'p Of^ i^^a-as^ (^7) 

reduces to the equation of Readhead (1994) for the equipartition Doppler factor (see eq. |A2|). 



The equipartition Doppler factor for a jet-like geometry is derived by using the 
appropriate expressions for this case, which are Ip = Sop (1 + z)^ ,5-2- [20/(5+20;)] ^-2^ 
Up = uop (1 + z) 5-i+[2/{5+2a)] ^ ^ ^ ^/O^a 9d aor{z) (1 + z)'^ 6'^ {sm(p)-\ where 9a and 
9b are angular sizes of the major and minor axes respectively, and <j) is the outflow angle (see 
GPCM93). The result using these substitutions is 

6eq{jetf^''' = Geqiaf GssAia)-' [aor{z)r' (1 + z)^''+'''^ S^'p 9f^ z^g^-^^) (sin,/.)^ (0„/0,), 

(A8) 

which reduces to equation (11). 

The brightness temperature, Tb, is a useful quantity to use when discussing radiative 
processes. The observed brightness temperature at the peak is 



So 

'd'^op 



Tbo = 1.77 X 10^2-2^ (A9) 



(e.g. Readhead 1994). The intrinsic brightness temperature, Tg,, can be related to the observed 
brightness temperature for a relativistic moving sphere using Tbi = Tbo{^ ~\~ z)/5, where z is the 
redshift to the source and 5 is the Doppler factor. These estimates are given by: 



TB^ = 1.77 X 1012 Sop (1 + z 



'd'^op 



(AlO) 



(e.g. Readhead 1994). Estimates of Tsi using 5eq are denoted TBi{eq), and those using djc are 
denoted TsiilG). 
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Table 1. Mean Brightness Temperatures 





Mean 


Mean 


Mean 


Class 


Tso/lO" K 


Tsi (eg)/10" K 


TBi(7C7)/10" K 


All Sources 


2.3 ±0.3 


0.87 ±0.03 


0.73 ± 0.03 


BL Lacs 


2.2 ±0.6 


0.86 ± 0.06 


0.75 ± 0.06 


CDHPQs 


3.5 ±0.6 


0.79 ± 0.03 


0.75 ± 0.06 


CDLPQs 


2.4 ±0.5 


0.89 ± 0.05 


0.80 ± 0.09 


CDQ-NPIs 


0.74 ± 0.24 


0.73 ±0.19 


0.51 ±0.13 


LDQs 


2.0 ± 1.7 


1.0 ±0.1 


0.78 ± 0.08 


LDQS-3C216 


0.30 ±0.09 


1.1 ± 0.1 


0.76 ± 0.08 


RGs 


0.36 ±0.10 


0.97 ±0.13 


0.65 ±0.13 
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Table 2. Mean 5eq, Sjc for True and Reassigned Angular Sizes 



Class No. True Reassigned 







Mean Seq 


Mean 6ic 


Mean 5eq 


Mean 6ic 


All Sources 


100 


6.0 ± 0.9 


5.1 ±0.6 


26 ±14 


9.2 ±2.4 


BL Lacs 


27 


4.5 ± 1.2 


3.4 ±0.7 


12 ± 7 


5.8 ±2.7 


CDHPQs 


24 


9.0 ± 1.5 


8.4 ±1.0 


90 ±56 


23± 10 


CDLPQs 


22 


6.9 ±1.5 


6.2 ± 1.4 


21 ± 12 


13 ±6 


CDQ-NPIs 


7 


5.1 ±2.2 


5.3 ±1.8 


26 ±18 


14 ±9 


LDQs 


11 


6.1 ±5.5 


3.6 ±2.9 


3.5 ±2.0 


2.0 ±0.9 


LDQS-3C216 


10 


0.56 ± 0.22 


0.66 ±0.22 


0.9 ±0.5 


0.6 ±0.2 


RGs 


9 


0.55 ±0.17 


0.65 ±0.14 


1.3 ±0.8 


1.1 ±0.5 
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Table 3. Mean Luminosity Densities and Luminosities 





Mean 


Mean 


Mean 


Mean 


Mean 


Mean 


Class 


log L^o 


logLo 


log L^i{eq) 


log Li (eg) 


logL.^/C) 


logL.(JC) 




( ergs \ 




( ergs \ 




/ ergs \ 








(^) 











All Sources 


33.6 ±0.1 


43.3 ± 0.1 


32.9 ±0.2 


42.7 ±0.3 


32.6 ±0.2 


42.2 ± 0.2 


BL Lacs 


32.8 ±0.2 


42.4 ± 0.2 


32.8 ± 0.5 


42.5 ±0.7 


32.6 ± 0.4 


42.3 ± 0.5 


CDHPQs 


34.3 ±0.1 


44.1 ±0.2 


32.0 ±0.3 


41.3 ± 0.4 


31.9 ±0.2 


41.1 ±0.3 


CDLPQs 


34.4 ± 0.2 


44.2 ± 0.3 


33.0 ±0.5 


42.6 ±0.6 


32.7 ±0.3 


42.2 ± 0.4 


CDQ-NPIs 


34.0 ± 0.3 


43.8 ± 0.4 


32.8 ± 0.7 


42.4 ±0.9 


32.5 ±0.5 


42.0 ±0.7 


LDQs 


33.2 ±0.3 


43.0 ±0.3 


34.4 ± 0.8 


44.8 ± 1.1 


34.0 ±0.7 


44.3 ± 0.9 


LDQS-3C216 


33.1 ±0.3 


42.9 ± 0.3 


35.0 ±0.7 


45.6 ±0.9 


34.5 ± 0.5 


45.0 ± 0.7 


RGs 


31.9 ±0.5 


41.7 ±0.5 


33.8 ± 1.1 


44.3 ± 1.3 


33.1 ±0.7 


43.3 ± 0.9 
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Table 4. The Overlap of the GPCM and VC Samples 



Source 


Name 




z 




Seq 


&IC 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


BL Lacs 














0454+844 







.112 


0.9 + 0.3 


1.3 


2.4 


0851+202 


OJ 287 


0, 


.306 


4.0 + 0.4 


11 


6.8 


1101+384 


Mkn 421 


0, 


.031 


2.50 + 0.04 


2.2 


0.92 


1308+326 







.996 


18 + 9 


4.0 


5.2 


1749+701 







.770 


10 + 1 


0.6 


0.9 


1803+784 







.684 


0.1 + 0.9 


5.3 


6.6 


2007+776 







.342 


3.4 + 0.7 


2.8 


3.6 


2200+420 


BL Lac 


0, 


.069 


4.6 + 0.2 


5.0 


3.4 


CDHPQ 














0212+735 




2, 


.370 


8 + 4 


5.1 


7.1 


0234+285 


CTD 20 


1, 


.213 


16 + 8 


6.6 


13 


1156+295 


4C 29.45 


0, 


.729 


41 + 2 


> 2.1 


> 4.9 


1253-055 


3C 279 


0, 


.538 


3.7 + 0.5 


13.2 


14.0 


1641+399 


3C 345 


0, 


.595 


9.6 + 0.2 


1.4 


4.1 


2223-052 


3C 446 


1, 


.404 


0.0 + 3.6 


16.6 


16.0 


2230+114 


CTA 102 


1, 


.037 


+ 24 


> 0.9 


> 1.5 


2251+158 


3C 454.4 





.859 


3.2 + 0.6 


> 5.2 


> 4.6 


CDLPQ 














0016+731 




1, 


.781 


16 + 4 


5.2 


7.9 


0153+744 




2, 


.340 


1.9 + 3.5 


> 1.3 


> 1.8 


0333+321 


NRAO 140 


1, 


.258 


8.3 + 0.5 


27 


13.0 


0430+052 


3C 120 


0, 


.033 


4.7 + 0.5 


> 11 


> 4.1 


0552+398 


DA 193 


2, 


.365 


3.6 + 1.8 


1.1 


2.2 


0711+356 


OI 318 


1, 


.620 


0.0 + 1.7 


20 


6.4 


0836+710 


4C 71.07 


2, 


.170 


16 + 3 


7.7 


6.7 


0923+392 


4C 39.25 


0, 


.699 


6.2 + 0.3 


6.5 


8.9 


1226+023 


3C 273 


0, 


.158 


8.8 + 0.2 


7.9 


4.6 


1928+738 


4C 73.18 





.302 


6.8 + 0.3 


3.2 


3.4 


CDQ-NPI 














0615+820 




0, 


.710 


1.7 + 1.7 


0.7 


1.4 


1039+811 




1, 


.260 


< 4 


16.6 


12.2 


1150+812 




1, 


.250 


6 + 3 


8.2 


9.6 
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Table 4 — Continued 



Source 


Name 




z 






Seq 


<5/c 


(1) 


(2) 


(3) 


(4) 






(5) 


(6) 


LDQ 


















0850+581 


4C 58.17 


1, 


.322 


7 + 


1 




2.2 


2.5 


0906+430 


30 216 





.670 


6± 


1 




61 


33 


1040+123 


3C 245 


1, 


.029 


5.2 ± 


2. 


,5 


0.5 


1.4 


1222+216 


40 21.35 


0, 


.435 


2.1 + 


0. 


9 


1.2 


1.0 


1618+177 


30 334 





.555 


2.8 ± 


0. 


,9 


> 0.2 


> 0.3 


1721+343 


40 34.47 


0, 


.206 


2.9 ± 


0. 


3 


0.1 


0.2 


1830+285 


40 28.45 


0, 


.594 


3.9 + 


0. 


8 


0.4 


0.4 


RG 


















0108+388 


GO 314 





.669 


0.90 + 


0. 


,25 


0.2 


0.7 


0316+413 


30 84 





.018 


0.6 + 


0. 


1 


1.2 


1.2 


0710+439 







.518 


0.1 ± 


0. 


2 


0.2 


0.4 


1228+127 


M 87 


0, 


.004 


0.2 ± 


0. 


1 


0.8 


0.8 


1637+826 


NGC 6251 


0, 


.023 


0.1 ± 


0. 


1 


> 1.3 


> 1.0 


2021+614 


OW 637 


0, 


.227 


0.3 + 


0. 


,3 


0.9 


1.1 


2352+495 


GZ 488 





.237 


< 


.4 




0.3 


0.5 
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Table 5. Solutions for T and (f) 



Source 


Name 


z 


Teg 


(deg) 




<^/c(deg) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


BL Lacs 














0454+844 




0.112 


1.4 + 0.2 


48 ± 97 


1.6 + 0.8 


18 + 35 


0851+202 


OJ 287 


0.306 


6 ± 8 


3 + 9 


5 + 2 


7± 11 


1101+384 


Mkn 421 


0.031 


3 ± 1 


26 ± 34 


5 + 4 


39 + 7 


1308+326 




0.996 


42 ± 72 


6 + 3 


33 + 40 


6 + 3 


1749+701 




0.770 


83 ± 133 


12 + 1 


56 ± 56 


12 + 1 


1803+784 




0.684 


3 + 4 


0.6 ± 4.5 


3 + 3 


0.4 ± 2.6 


2007+776 




0.342 


3.6 ± 1.6 


21 ± 26 


3.5 ± 0.7 


16 ± 17 


2200+420 


BL Lac 


0.069 


4.7 ± 0.4 


12 ± 20 


5.0 ± 1.6 


16 ± 11 


CDHPQ 














0212+735 




2.370 


9 + 9 


10 ± 10 


8 + 5 


8 + 7 


0234+285 


CTD 20 


1.213 


23 + 32 


6 + 4 


16 + 10 


4 + 3 


1156+295 


4C 29.45 


0.729 


< 405 


< 2.8 


< 175 


< 2.8 


1253-055 


3C 279 


0.538 


7 + 10 


2 + 7 


7.5 + 6.5 


2 + 4 


1641+399 


3C 345 


0.595 


34 + 52 


12 ± 1 


13 + 9 


10 + 3 


2223-052 


3C 446 


1.404 


8 + 13 


0.0 + 1.5 


8 + 8 


0.0 + 1.6 


2230+114 


CTA 102 


1.037 


> 1.0 


+ 180 


> 1.1 


+ 180 


2251+158 


3C 454.4 


0.859 


> 3.7 


< 10 


> 3.5 


< 12 


CDLPQ 














0016+731 




1.781 


27 ± 36 


7 + 2 


20 + 14 


6 + 2 


0153+744 




2.340 


> 2.1 


< 42 


> 2.1 


< 33 


0333+321 


NRAO 140 


1.258 


15 + 20 


1 + 4 


9 + 4 


4 + 6 


0430+052 


3C 120 


0.033 


> 6 


< 4 


> 4.5 


< 14 


0552+398 


DA 193 


2.365 


7 + 11 


29 + 14 


4.2 + 3.5 


23 + 13 


0711+356 


OI 318 


1.620 


10 + 16 


+ 1 


3 + 3 


+ 5 


0836+710 


4C 71.07 


2.170 


19 + 20 


6 + 4 


21 ± 16 


6 + 2 


0923+392 


4C 39.25 


0.699 


6.3 + 0.5 


9 + 15 


6.6 + 2.3 


6 + 8 


1226+023 


3C 273 


0.158 


8.9 + 1.5 


7+10 


11+6 


10 + 4 


1928+738 


4C 73.18 


0.302 


8.9 + 9.1 


14 + 8 


8.6 + 5.2 


14 + 5 


CDQ-NPI 














0615+820 




0.710 


3 + 6 


54 + 42 


2 + 2 


41 ± 31 


1039+811 




1.260 


< 8.8 


< 1.5 


< 6.7 


< 2.7 


1150+812 




1.250 


6 + 4 


7+14 


7 + 3 


5 + 8 
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Table 5 — Continued 



Source 


Name 




z 


Teg 


<*<!<7 ((leg) 




0/c(deg) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


LDQ 
















0850+581 


4C 58.17 


1, 


,322 


12 ± 16 


15 + 5 


11 ± 10 


15 + 4 


0906+430 


3C 216 


0, 


.670 


31+49 


0.2 + 0.6 


17+16 


0.6 + 1.2 


1040+123 


30 245 


1, 


.029 


28 + 51 


22 + 10 


11 + 13 


21 + 9 


1222+216 


4C 21.35 


0. 


.435 


2.8 + 2.9 


41 + 30 


3.2 + 2.9 


44 + 18 


1618+177 


3C 334 


0, 


.555 


> 3 


< 39 


> 3 


< 39 


1721+343 


4C 34.47 


0, 


,206 


46 + 74 


38 + 4 


32 ± 32 


38 + 4 


1830+285 


4C 28.45 


0, 


,594 


20 + 32 


29 + 6 


21 ± 23 


29 + 6 


RG 
















0108+388 


OC 314 


0, 


,669 


4 + 6 


94+17 


1.7+1.0 


81 + 32 


0316+413 


30 84 


0, 


,018 


1.15 + 0.07 


60 + 160 


1.15 + 0.08 


56 + 99 


0710+439 




0, 


,518 


2.6 + 3.9 


170 + 25 


1.5 + 1.1 


170 ± 30 


1228+127 


M 87 


0, 


,004 


1.0 + 0.4 


120 + 180 


1.1 ± 0.3 


130 ± 90 


1637+826 


NGC 6251 


0, 


,023 


> 1.1 


< 15 


> 1.0 


< 88 


2021+614 


OW 637 


0, 


,227 


1.1 + 0.3 


110 + 180 


1.04 + 0.09 


60 + 180 


2352+495 


OZ 488 


0, 


,237 


< 1.8 


> 133 


< 1.25 


> 124 
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Table 6. Median Values of T and (p 



Class 


Med. Feq 


Med. (peqideg) 


Med. Tic 


Med. (picideg) 


BL Lacs 


4.2 ±0.7 


12 ±4 


5.0 ±0.8 


14 ±2 


CDHPQ 


9± 1 


4±2 


8±1 


3.0 ±1.5 


CDLPQ 


9.5 ± 1.2 


7±2 


8.8 ± 1.7 


6± 1 


CDQ-NPI 


4.5 ± 1.5 


30 ±23 


4.5 ±2.5 


23± 18 


LDQ 


24 ±5 


26 ±6 


14 ±3 


25 ±6 


RG 


1.15 ±0.08 


110 ±9 


1.15 ±0.06 


81 ±14 
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Fig. 1. — (a) The equipartition Doppler factor, Seq, vs. (1 + z); (b) the inverse Compton Doppler 

factor, 5ic, vs. (1 + z) for the 100 sources in this study. Symbols: sohd circles are BL Lacs, solid 
diamonds are CDHPQs, solid squares are CDLPQs, solid triangles are CDQ-NPIs, open diamonds 
are LDQs, and open squares are RGs. 

Fig. 2. — The distribution of equipartition Doppler factors, Sgq (dotted line), and inverse Compton 
Doppler factors, Sjc (dashed line), for all sources in this study. 

Fig. 3. — The distribution of equipartition Doppler factors, 5eq (dotted line), and inverse Compton 
Doppler factors, 5jc (dashed line), for each class of AGN: (a) BL Lacs, (b) CDHPQs, (c) CDLPQs, 
(d) CDQ-NPIs, (e) LDQs, (f) RGs. 

Fig. 4. — The distribution of observed brightness temperatures, Tbo (solid line), and intrinsic 
brightness temperatures, TBi{eq) (dotted line) &; TBi{IC) (dashed line), for all sources in this 
study. 

Fig. 5. — The distribution of observed brightness temperatures, Tbo (solid line), and intrinsic 
brightness temperatures, TBi{eq) (dotted line) & TBi{IC) (dashed line), for each class of AGN: (a) 
BL Lacs, (b) CDHPQs, (c) CDLPQs, (d) CDQ-NPIs, (e) LDQs, (f) RGs. 

Fig. 6. — The uncorrected luminosity density, L^^^, vs. (1 -|- z) for all 100 sources. Symbols: solid 
circles are BL Lacs, solid diamonds are CDHPQs, solid squares are CDLPQs, solid triangles are 
CDQ-NPIs, open diamonds are LDQs, and open squares are RGs. 

Fig. 7. — The uncorrected luminosity , vs. (l+z) for all 100 sources. Symbols: solid circles are 
BL Lacs, solid diamonds are CDHPQs, solid squares are CDLPQs, solid triangles are CDQ-NPIs, 
open diamonds are LDQs, and open squares are RGs. 

Fig. 8. — The intrinsic luminosity density Li,i{eq) vs. the rest frame frequency Vi{eq). The result of 
a power law fit is Lui{eq) oc Ui{eq)^'^^^'^ with a reduced equal to 4.6 and a correlation coefficient 
r = 0.91. Symbols: solid circles arc BL Lacs, solid diamonds arc CDHPQs, solid squares are 
CDLPQs, solid triangles are CDQ-NPIs, open diamonds are LDQs, and open squares are RGs. 

Fig. 9. — The intrinsic luminosity Li{eq) vs. the rest frame frequency Ui{eq). The result of a power 
law fit is Li{eq) oc Ui(eqY''^^^'^ with a reduced equal to 4.0 and a correlation coefficient r = 0.96. 
Symbols: solid circles are BL Lacs, solid diamonds are CDHPQs, solid squares are CDLPQs, solid 
triangles are CDQ-NPIs, open diamonds are LDQs, and open squares are RGs. 

Fig. 10. — The intrinsic luminosity density Lyi{eq) vs. (l + z). The result of a power law fit 
is L[^i{eq) oc (1 -H z)'^-^^^-^ with a reduced equal to 1.3 and a correlation coefficient r = 0.16. 
Symbols: solid circles are BL Lacs, solid diamonds are CDHPQs, solid squares are CDLPQs, solid 
triangles are CDQ-NPIs, open diamonds are LDQs, and open squares are RGs. 

Fig. 11. — The intrinsic luminosity Li{eq) vs. {l + z). The result of a power law fit is 
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L,yi{eq) oc (1 + 2)0-9±2.2 ^j^j^ reduced equal to 1.4 and a correlation coefficient r = 0.04. 
Symbols: solid circles arc BL Lacs, solid diamonds are CDHPQs, solid squares are CDLPQs, solid 
triangles are CDQ-NPIs, open diamonds are LDQs, and open squares are RGs. 

Fig. 12. — The apparent speed in the plane of the sky, Papp, vs. redshift, z for the overlap of the 
GPCM93 and VC94 samples (upper bounds are indicated by arrows). Symbols: solid circles are 
BL Lacs, solid diamonds are CDHPQs, solid squares are CDLPQs, solid triangles are CDQ-NPIs, 
open diamonds are LDQs, and open squares are RGs. 

Fig. 13. — Teq vs. (f>gq, estimates using the equipartition Doppler factor for the overlap of the 
GPCM93 and VC94 samples (upper and lower bounds indicated by arrows): (a) with errors shown, 
and (b) without errors shown. Symbols: solid circles are BL Lacs, solid diamonds are CDHPQs, 
solid squares are CDLPQs, solid triangles are CDQ-NPIs, open diamonds are LDQs, and open 
squares are RGs. 

Fig. 14. — Fjc vs. (pic, estimates using the inverse Compton Doppler factor for the overlap of the 
GPCM93 and VC94 samples (upper and lower bounds indicated by arrows): (a) with errors shown, 
and (b) without errors shown. 

Fig. 15. — (a) (peg and (b) Tgg as functions of {1 + z). Symbols: solid circles are BL Lacs, solid 
diamonds are CDHPQs, solid squares are CDLPQs, solid triangles are CDQ-NPIs, open diamonds 
are LDQs, and open squares are RGs. 

Fig. 16. — (a) (pic and (b) Tic as functions of {1 + z). Symbols: solid circles are BL Lacs, solid 
diamonds are CDHPQs, solid squares are CDLPQs, solid triangles are CDQ-NPIs, open diamonds 
are LDQs, and open squares are RGs. 
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